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FOREWORD 


The  two  objectives  of  Contract  A.F  29(601)-4283,  (SRI  Project 
No.  PU-3514)  are: 

(1)  The  devolopraont  of  uncouplocl  wave  equations  for  a  nonhomoge- 
neous  nxisymmotrlc  half-space,  suitable  for  treatment  by  the 
method  of  characteristics. 

(2)  The  theoretical  study  using  the  method  of  characteristics  of 
stress,  strains,  and  particle  motions  in  layered  and  non- 
homogeneous  elastic  half-spaces  which  are  surface  loaded  by 
a  nuclear  airblast. 

This  contract,  from  January  1961  to  June  1962,  was  a  continuation  of  a 
previous  study  (Contract  AF  29(601)-19'18 ,  SRI  Project  No.  SU-2893) . 

The  research  effort  is  divided  into  five  integral  phases  which  will 
be  reported  separately.  This  is  the  final  report  for  Phase  I  which  is 
the  development  of  uncoupled  wave  equations  for  a  nonhomogeneous  axi- 
symmetric  half-space.  Phase  II  is  the  development  of  the  numerical 
calculation  method  for  a  homogeneous  and  a  layered  half -space.  Phase  II 
includes  coding  and  check-out  of  the  resulting  program.  Tho  final 
report  on  this  phase  will  pi'osent  detailed  explanations  of  mathematical 
mothods  developed  and  sample  calculations.  Production  runs  and 
analysis  of  rosults  for  u  layered  half-space  will  constitute  Phase  III. 
This  phase  will  establish  tho  influence  on  ground  motion  of  variations 
in  independent  problem  parameters,  o.g.,  tlio  influence  on  the  maximum 
stress  field  duo  to  the  difference  in  comprossional  wave  speeds  between 
layei’s,  to  depth  of  overburden,  and  to  dlfforenoos  in  Poisson's  ratio. 

A  complete  description  of  tho  codes  used,  along  with  Instructions  to  the 
user,  will  bo  presented  at  that  time.  Production  runs  and  analysis  of 
results  for  a  nonhomogeneous  half-space  will  constitute  Phase  IV.  In 
Phase  V  theoretical  I’os.ults  will  bo  compared  with  field  data.  In  this 
phase,  tho  rosults  of  Phases  III  and  IV  will  bo  compared  with  available 
tost  data  and  final  conclusions  pi’osuntod  as  to  tho  applicability  and 
limitations  of  tho  theoretical  results. 
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ABSTHACT 


Uncoupled  wave  equations  are  developed  lor  a  nonhomogeneous  half- 
space  for  two  cases,  (1)  where  Poisson's  ratio  is  constant,  and  (2)  where 
Poisson's  ratio  increases  with  depth  in  a  nonarbitrary  but  realistic 
manner.  The  resulting  wave  equations  are  numerically  solvable  by  the 
method  of  %'jharacteri8tic8  since  they  are  identical  to  those  for  a 
homogeneous  medium. 
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1. 


INTRODUCTION 


Tills  report  describes  the  work  toward  obtaining  separable  equa¬ 
tions  dT  motion  which  govern  the  propagation  of  elastic  disturbances 
in  nonhomogenoous  media.  Separability  here  means  two  Independent 
equations  of  motion  which  may  be  subject  to  related  boundary  conditions. 
Attention  has  been  confined  to  the  case  of  axial  symmetry  where  the  Larn^ 
parameters  (^, P )  and  the  donslty  (p)  are  allowed  to  vary  only  as  func¬ 
tions  of  the  distance  below  the  surface  of  the  half-space. 

1  + 

This  problem  was  first  analyzod  by  Hook  who  was  able  to  treat  five 
special  cases  (Table  I),  all  of  which  require  a  specific  relationship 
between  p,  p,  and  Poisson's  ratio,  a.  Hook's  approach,  however,  suggests 
the  possihillLy  of  a  more  general  approach  wherein  such  restrictions  are 
not  required.  The  practical  Importance  of  obtaining  a  generalized  sep¬ 
aration  of  the  equations  of  motion  Is  that  this  would  allow  direct 
nppllcatlou  of  the  calculation  prt>cedui'es  developed  in  Part  I  of  this 
report  to  cases  whore  there  Is  a  continuous  (as  i -'posed  to  a  stepwise) 
variation  of  elastic  properties  with  depth. 


TAUU3  I 

SUMMARY  OF  NONIIOMOfiRNKOlIS  CASKS  TllKATEU  IN  REF,  NO.  1 


CASE 

Vai-iatlon  with  depth,  z  | 

Shear  Modulus 

Density 

Poisson's  Ratio 

1 

p  constant 

-2 

/>  ■/, 

(r  constant 

2 

n 

n-2 

p  <x  z 

O’  =  1/(2  +  n) 

■I 

i\ 

fX  tX  Vi 

li(u-l) 

/)  <•'  z 

cr  -  (3  -  2n)/4 

H 

n 

H  7. 

n-1 

p  It  z 

(T  =  l/(n  -i  1) 

IB 

-  1 
fx 

p  arbitrary 

t;  constant 

Note  that  I’olssou's  I'utlo  is  an  arbitrary  constant  (Cases  1  and  5 
la  Table  I)  or  a  constant  wlilch  depends  on  tile  variation  of  shear  modu¬ 
lus  with  dojith,  l.e,,  which  depends  on  n.  Hence,  the  cumpresslunal  Wave 
velocity.  O',  varies  us  the  shear  wave  velocity,  p.  Tho  loss  in  gener¬ 
ality  results  from  tho  unusual  variation  ol  either  p  or  p  (Cases  1 
and  5);  (.ir  if  p  is  constant,  then  (t  ==  l/i  (Cases  2  and  3)  or  a  =  1/2 
(Case  4).  For  (T  =  1/2,  the  shear  eontrlbutlon  may  bo  nogloeted  Initially, 
In  which  case  the  variation  of  compress lonal  wave  velocity  with  depth 
may  be  chosen  arbitrarily.  This  loaves  Cases  2  and  3  of  practical  im¬ 
portance,  with  />  t'oiiMiant  and  (r  =  l/'l  (Case  2,  n  ~  2;  Case  3,  n  -  1). 

The  main  mathomuticui  result  of  this  investigation  is  that  sopara- 
blltty  may  be  obtained  provided  tliai  A  ,  p  ,  and  p  are  related  by  tho 
ordinary  d  llT'oi-ont  lai  equations 


+  Numbered  references  are  listed  at  the  end  of  this  report. 
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(A  +  2H)  d^fi  _  /dju\ 


(A  I  2fi)  dp 

2p  dz  ~  dz 

Solutions  ol'  those  equations  loud  to  physically  acceptable  specifica¬ 
tions  for  tho  variation  of  elastic  parameters  with  depth. 

The  results  of  tho  Investigation  are  examined  In  more  detail  in 
Section  2,  while  a  method  by  which  these  results  are  obtained  is  given 
in  Section  3. 


2.  RESULTS  AND  EXAMPIJ<;S 

Tho  differential  equations  which  govern  the;  propagation  of  elastic 
waves  in  nonhomogencous  as  well  as  homogoneous  media  are  the  two  equa¬ 
tions  of  dynamic  equilibrium  together  with  tlie  four  equations  given  by 
the  generalized  Hooko's  law 
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fbu  u\ 

(or  r) 

(5) 

(T 

rz 

/Dll 
=  ^{d2 

^  Sr^ 

(6) 

In  thuso  equal  ions  tr  ^  ,  (X^ ,  tr^,  ,  u,  and  w  ui’o  tlio  usual  eomponents  of 
stress  and  displacement  when  eonditleiis  of  axial  symmetry  prevail,* 
and  A  is  tlie  dilutution  given  by  ^  i 


*  See,  for  example,  Timoshenko,  S. ,  Theory  of  Klustlcity,  MeGi’aw  Hill 
U<jok  Co.,  Inc.,  Now  York  and  London,  193'1. 


It  is  instructive  to  consider  some  examples  In  order  to  investigate 
to  what  extent  the  above  results  may  bo  applicable  in  some  possible  and 
practical  cases. 

Example  1:  Constant  Poisson's  Ratio 
We  define 


y 


(X  +  2n) 


Thus,  in  terms  of  Poisson's  ratio, 


y 


2(1  -  Cl) 

(1  -  2ff) 


and 


cr 


1  C2  -  y) 

2  (1  -  y)  • 


It  is  evident  that  if  0  £  a  £  then  2  <  y  s  oo  .  In  terma  of  y  It  Is 
convenient  to  rewrite  Bcjuation  (12)  in  the  form 


We  consider  first  the  case  cr  =  constant.  In  this  case  Equation  (12) 
is  readily  Integrated  to  yield 


H  = 

whore  A  and  B  are  arbitrary  constants.  The  density,  p,  is  then  found  by 
integrating  Equation  (11)  or  (1.3)  to  give 


(t^,)  " 


1  (l-o')/cr 


z  +  A 


(15) 


p  -  C 


(r^)  ' ' »] 


(l-2ff)/(T 


(16) 


where  C  is  an  arbitrary  constant.  Since  any  numerical  Integrals  of  the 
equations  of  motion  can  at  best  bo  carried  out  for  only  a  finite  time 
and  hence  a  finite  range  of  z,  Equations  (15)  and  (1.6)  can  probably  be 
made  to  fit  many  practical  cases  whore  only  a  finite  depth  of  the  medium 
is  of  interest. 


Equations  (15)  and  (16)  may  bo  put  Into  a  mure  convenient  form  for 
computation.  Define  p(0)  and  p(0)  as  tho  density  and  shear  modulus  at 
zero  depth;  then 


(l-cr)/(T 


jfe)  -  "  '1 

(15a) 

(1-2(t)/(7 

jfe  ■  h«i'>  * '] 

(16a) 
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where  and  C2  are  new  arbitrary  constants.  From  Equations  (15a)  and 
(16a)  it  may  be  seen  that  the  density  and  shear  modulus  are  power  law 
variations  with  depth  and  that  Cj  is  a  scale  factor  which  adjusts  the 
depth  at  which  the  shear  modulus  attains  a  given  value.  Typical  vari¬ 
ations  with  depth  of  density,  shear  modulus,  and  shear  wave  velocity 
which  may  be  obtained  using  Equations  (15a)  and  (16a)  are  illustrated 
in  Fig.  1. 


Example  2;  Variable  Poisson's  Ratio 


We  consider  next  the  case  whoi’o  Poisson's  ratio  varies  from  some 
arbitrary  positive  value  less  than  1/2  to  1/2  as  0  —  z  — ‘®.  In  par¬ 
ticular  we  assume 


a 


(z  +  b)  -a 

- T~S 

2(z  +  b)  -a 


or 


y 


2(z  +  b)’^ 

2 

a 


(17) 


The  first  integral  of  (12)  gives 


iL 

M’ 


(z  +  b) 


+  C, 


where  is  an  arbitrary  constant  of  integration.  Then  with  2a  +  b 


+  a  +  b)oxp 


a 

z  +  a  +  b 


(18) 


where  Cg  is  an  arbitrary  constant.  The  density,  p,  can  then  bo  found 
from  (13)  with  the  aid  of  (18)  to  give 


P 


r,  (2  +  b) 

“3  (z  +  a  +  b) 


exp 


_ a _ 

(z  +  a  +  b) 


(19) 


where  C3  Is  an  arblti'ary  constant.  It  should  be  noted  that  p  as  defined 
by  (18)  increases  without  limit,  while  the  density  increases  to  an 
asymptotic  constant  value  with  increasing  depth.  The  variations  of  the 
speeds  of  the  so-called  dllatatioiml  and  rotational  waves  with  depth  are 
easily  written  down  from  (17),  (18)  and  (19)  to  give 


and 


2 


a 


2  —  (z  +  a  +  b) 
^3 


2(z  +  b) 

2 

a 


^2  (z  +  a  +  b)^ 
Cg  (z  +  b)' 


(20) 


(21) 
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DENSITY,  />,  AND  SHEAR  WAVE  VELOCITY,  0 


FIG.  1  TYPICAL  DLNSITY,  SHEAR  MODULUS,  AND  SHEAR  WAVE  VELOCITY  PROFILES 
FOR  CONSTANT  POISSON’S  RATIO, 

Cl  «  1/4  [;i(0),  p(0),  and  /i(0)  are  values  of  parometers  at  lero  depth] 
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Again,  these  variations  can  be  made  to  lit  many  possible  cases  for  finite 
depths.  Equation.s  (17),  (19),  (20)  and  (21)  may  be  written  in  the 
following  forms: 


(1  +  b  z)^  -  a  ^ 

- 5:  (17a) 

2(1  +  b^z)^  - 


a 


Q’(O) 


(1  +  aj^)(l  +  b^z) 


oxp 


1  +  a^  +  bj^z 


(1  +  +  b,z) 

- ^ - i—  Vl  H-b^Z 

1  +  a,  1 


(1  +  +  bj^z)  1 


^  +  '‘l  Vl  +  b^z 


(19a) 


(20a) 


(21a) 


whore  b^  Is  an  arbitrary  constant  serving  as  a  depth  scaling  factor  and 
is  determined  by  the  Poisson's  ratio  at  zero  depth,  i.e., 


or 


0(0) 


(17b) 


1  -  2p 
1  -  P 


(17c) 


Typical  variations  with  depth  obtained  using  Equations  (17a),  (19a), 
(20a),  and  (21a)  are  illustrated  in  Pig.  2  for  (7(0)  =  0.265. 

The  above  results  do  not  exhaust  all  possible  cases  for  variable  (7; 
however,  each  desired  variations  of  a  would  have  to  be  substituted  Into 
Equations  (11)  and  (12)  and  investigated  separately. 
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FIG.  2  TYPICAL  DENSITY,  POISSON'S  RATIO,  COMPRESSION  WAVE  VELOCITY, 

AND  SHEAR  WAVE  VELOCITY  PROFILES  FOR  VARIABLE  POISSON'S  RATIO 

[MO),  01(0),  and  /i(0)  are  values  of  parameters  at  zero  depth] 
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3.  SEP'ARATION  OF  TJIE  EXiUATIONS  OF  MOTION 

When  Equations  (3)  to  (6)  aro  used  to  eliminate  the  stresses  in 
Equations  (1)  and  (2),  the  resulting  equations  which  take  into  account 
the  z  dependence  of  X,  ,  and  p  are 


_d_ 

8r 


|(X  +  2p)[g^  +  7  5^ Ir  aife  "  8?)|  ^'^  8?  =  P  -1 

ot 


and 


_8_ 

8z 


f/\  ,  o  u  8w\l  1  8  /8u  8w\  „8  /8u  u\  8^w 

[(X  +  2p)(_  .  -  +  g-)]  +  -  ^rp(-  +  -  2^p(^  "  r)  = 


In  a  manner  similar  to  that  employed  In  Reference  1,  the  following 
substitutions  are  now  made  in  (17)  and  (18); 


8^ 

8r 


(f  *  »•) 
•  "  (S'  *  '*)  •  7 


(24) 


where  g  and  f  are  functions  of  z  alone.  Equations  (22)  and  (23)  then 
become 


8r 
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i  (pB  -  p')  -^1  =  0 

8t^ 


(22a) 


and 
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^  |(A  +  +  [(X  +  +  2/J']  ^  +  [(X  +  fi)f'  +  -  P  •— j 

+  (US  -  -  2p"  H  +  [pf"  -  +  (P  '  -  PS)  ~ 

^  dt^ 

r 

r 

+  [pg'  +  2n'g  -  ((X  +  2p)K)*]0  -P^j  =  0 

8t  ' 


(23a) 


Soparabillty  will  have  boon  aohioved  If  f  and  g  can  be  chosen  so 
that  Equations  (22a)  and  (23a)  can  be  put  into  the  form 


flr 


S  "  -  T.\  S 


—  +  P(z)j  Kj  +  7  -  0 


(25) 


who  re  Kl  and  K2  are  dii'foi’ontlal  expressions  involving  <t>  and  6  alone, 
respectively.  If  (22a)  and  (23a)  are  reduced  to  (25)  then  the  solution 
=  0  and  K<j  =  0  yields  two  separate  partial  differential  equations 
for  (p  and  0  , 


One  solution,  and  to  date  the  only  one  which  lins  been  obtained,  is 
f  =  0,  g  =,-^  ,  where  X,  p,  and  p  satisfy  (11)  and  (12).  The  separatetl 
equations  a/o  then  given  by  (9)  and  (10). 


d.  UISCUSSION 

The  work  reported  here,  and  especially  the  special  case  for  variable 
Poisson's  ratio  given  in  Section  2,  is  capable  of  representing  a  physical 
situation  which  Is  very  likely  .similar  to  the  actual  behavior  of  the 
clastic  constants  with  depth  in  the  earth.  Per  instance,  it  i.s  very 
likely  that  the  earth  becomes  more  incompressible  with  depth;  this  wovilU 
indicate  that  Poisson's  ratio  increases  with  depth  from  some  value  loss 
than  1/2  to  1/2  as  it  does  in  the  exuiiipie.  tiirthor,  it  is  unlikely  that 
the  mass  density  inewja.sus  withoiit  bound  with  depth;  the  real  behavior 
is  very  probably  a  gradual  increase  with  depth  to  an  asymptotic  constant 
value  as  given  by  Equation  (19).  The  resistance  of  the  earth  to  volume 
change  increases  more  rapidly  than  rosistaneo  to  shear  as  indicated  by 
the  resulting  expressions  foi- (V^  and In  Soctloii  2.  The  situation 
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where  Poisson's  ratio  is  constant  is  a  special  case  in  the  present  work, 
but  again  it  yields  a  space  variation  of  the  density  which  is  much  slower 
than  the  variation  of  either  of  the  wave  speeds. 

The  present  work  can  be  i-eadily  applied  to  the  general  three- 
dimensional  case  having  medium  property  variation  in  one  dimension  with 
only  sliglit  modifications.  It  is  the  opinion  of  the  authors  that  a  more 
general  solution  can  be  found  using  the  present  methods.  This  more 
general  solution  should  contain  all  of  the  results  in  Reference  2  as 
special  cases. 

Another  investigation  which  has  not  been  completed  but  which  shows 
promise  is  directed  toward  obtaining  more  general  differential  equations 
of  motion  than  those  reported  here  or  in  Reference  1.  It  Is  to  be  noted 
that  the  differential  equations  (9)  and  (10)  closely  resemble  those  which 
govern  the  propagation  of  waves  in  a  homogeneous  medium.  The  more  general 
equations  which  might  bo  obtained  are  of  the  form 

(p  (p 

Ip  +  a  “  +  b(p  +  cd)  +  dd>  -  — ^  =  0  (26) 

^rr  r  ^  2 

c 

where  b  and  are  such  that  the  characteristic  equation 
c 


2  2 
f  +  bf 
r  z 


-if  " 

2  t 


=  0 


(27) 


admits  real  solutions.  This  type  of  general Izat ion  should  permit  more 
freedom  in  the  space  variation  of  A ,  p. ,  and  p . 
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